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The Vis polyoxovanadate molecule is made of 15 spins 
1/2 with antiferromagnetic couplings. It belongs to the class 
of molecules with very large Hilbert space dimension (2 in 
V15, 10 8 in Mni2-ac). It is a low spin/large molecule with spin 
S = 1/2. Contrary large spins/large molecules of the Mni2-ac 
type, Vis has no energy barrier against spin rotation. Magne- 
tization measurements have been performed and despite the 
absence of a barrier, magnetic hysteresis is observed over a 
timescale of several seconds. This new phenomenon charac- 
terized by a "butterfly" hysteresis loop is due to the effect of 
the environment on the quantum rotation of the entangled 15 
spins of the molecule, in which the phonon density of states 
is not at its equilibrium (phonon bottleneck). 
PACS numbers: 75.50.Xx, 75.45. +j, 71.70.-d 
keywords: magnetic quantum effects, relaxation phenomena, 
dissipative two-level system. 



1. Introduction 

The recent evidence for quantum tunneling of the mag- 
netization in big molecules with large spin S = 10 M 
allows a link to be made between mesoscopic physics 
and magnetism In these systems, energy barriers 
are high and tunnel splittings between symmetric and 
antisymmetric quantum states are exponentially small 
(1CP 7 — 1CP 11 K). As a consequence, quantum relaxation 
is slow and only related to the spin baths |3|-|j| • Resonant 
phonon transitions are irrelevant, unless between states 
at different energies |(| or in the presence of a transverse 
field large enough to create a tunnel splitting of the or- 
der of the temperature energy scale Q. Here we show 
that the reversal of an ensemble of non-interacting spins 
S = 1/2 is irreversible as a consequence of different cou- 
plings to the environment. Such an effect gives rise in our 
system to a "butterfly" hysteresis loop. Spin rotation in 
V15 may be viewed as an experimental realization of the 
theoretical problem of a 2-level system with dissipation 
|^|,^) . Spin-phonon transitions take place between the two 
splitted ground states at well defined energies fko — Ajj, 
where Ah is the field dependent energy separation of the 
two levels. The width Aw of such transitions being rela- 
tively narrow, the number of available phonon states at 
energy fiu is much smaller than the number of molecules 
and a hole must be burned in the spin-phonon density of 
states jnj[ll). This hole, moving with the sweeping field, 
leads to a non-equilibrium spin-phonon density of states, 
giving in the V15 complex the observed phonon bottle- 



neck and hysteresis . A phonon "bottleneck plateau" , 
identified on the measured hysteresis loop, gives a char- 
acteristic butterfly shape. These results are corroborated 
with numerical calculations based on this model. 

2. The V15 magnetic moment evolution in low 
sweeping rates. 

The V15 complex of formula K e [V^AseO i2 (H 2 0)} ■ 

It is made 



8H2O was prepared as described in 12 
of a lattice of molecules with fifteen V 



ions of spin 

S = 1/2, placed in a quasi-spherical layered struc- 
ture formed of a triangle, sandwiched by two hexagons 
(Fig. Q). All the interactions being antiferromagnetic, 
the resultant (collective spin) is S = 1/2. The symmetry 
is trigonal (space group i?3c, a = 14.029 A, a = 79.26°, 
V = 2632 A 3 ). The unit-cell contains two V15 clusters 
and it is large enough so that dipolar interactions be- 
tween spin 1/2 molecules are negligible (few mK). Each 
hexagon contains three pairs of strongly coupled spins 
(J = —800 K) and each spin at a corner of the inner 
triangle is coupled to two of those pairs (one belonging 
to the upper hexagon and one belonging to the lower 
hexagon (J' J x -150 K, J" = J 2 = -300 K @). 
Thus, the V15 molecule can be seen as formed by three 
groups of five V IV ions with resultant spin 1/2 and as- 
sembled on each corner of the inner triangle. These three 
spins 1 /2 interact to each others through two main paths, 
one passing by the upper hexagon and one passing by 
the lower one. This is a typical example of frustrated 
molecule, where the exchange Jo between the spins 1/2, 
is much smaller than the exchange interactions between 
two spins. 

Two types of single crystal magnetization measure- 
ments have been performed: (i) characterization of the 
general thermodynamical properties of the system at 
equilibrium and (ii) study of the non-equilibrium to equi- 
librium transition. 

In the first case a small dilution refrigerator allowing 
measurements above 0.1 K, was inserted in an extrac- 
tion magnetometer providing fields up to 16 T, with low 
sweeping rates. Below 0.9 K we observed one jump in 
zero field and two others at B x = ± 2.8 T (Fig. |). 
They correspond to the |l/2, -1/2) |l/2, 1/2) and 
|l/2, 1/2) <-> |3/2,3/2) spin transitions, with respective 
saturations at 0.50 ± 0.02 fiB and 2.95 ± 0.02 hb per 
V15 molecule, in agreement with the level scheme given 
in Fig. |^ inset. For all temperatures the magnetization 
curves are reversible and do not show any anisotropy, 



1 



showing that eventual energy barrier preventing spin re- 
versal must be quite small (less than 50 mK). 

The magnetization curves measured at equilibrium 
(Fig. H|) are fitted using the Heisenberg Hamiltonian for 
three frustrated spins 5 = 1/2: 



H- — Ja(S a lS a 2 + S a 2S a 3 + S a 3S a i) 

-gH B Bo($i+& + $a) (0.1) 



a=X,Y,Z 



with S l = S 2 = S 3 = 1/2, J x = Jy < 0, J z < (antifer- 
romagnetic coupling), g = 2, B the applied field and p,B 
the Bohr magneton. In the isotropic case (Jx,y,z = Jo) 
the configuration S — 3/2 becomes favourable against 
the S — 1/2 one for B\ = — 3 Jo/(2(?/xb). The measured 
value B\ = 2.8 T gives Jo = —2.5 K. The comparizon 
of the calculated curves (Szi — (i\Szi + Sz2 + Szz\i) 
averaged over the different eigenstates Ei,\i),i = 1..8) 
with the measured ones is quite good at any tempera- 
ture. The width of the transition at B\ = ±2.8 T is 
nevertheless broader in the experiments (Fig. ^[). This 
broadening of about 0.7 T cannot be inferred to dipo- 
lar or hyperfine field distributions (about 1 mT and 
40 mT respectively). The observed anisotropy of the 
electronic g-factor (<7 a =<ft,=1.95 and g c =1.98) is also too 
small |l4|] . Antisymmetrical Dzyaloshinsky-Moriya inter- 
actions ( jl5 16 1 and ref. therein) allowed by symmetry 
and coupling the states 5=1/2 and S = 3/2, could 
explain this broadening. 

In low fields, susceptibility measurements give the 
following effective paramagnetic moments: A*e// = 
g v / S(S+ 1) = 1.75 ± 0.02 /i B corresponding to S = 1/2 
below 0.5 K, and £t e // = 3 ± 0.02 /zb corresponding to 
3 independent spins 5=1/2 above 100 K. It confirms 
clearly that the ground state of the molecule is 5 — 1/2 
below 0.5 K. 

3. Spin-phonon coupling in V15 in low fields 
and high sweeping rates: a dissipative spin 1/2 
two-level system. 

We showed that, between -2.8 T and 2.8 T, the to- 
tal spin of this molecule is 5 = 1/2. Such a small spin 
leads to vanishingly small energy barrier and relatively 
large splitting in zero-field (~ 10~ 2 K). Resonant spin- 
phonon transitions between the symmetric and antisym- 
metric states at resonance are then possible in sub-Kelvin 
experiments, contrary to the case of big molecules in zero 
field, where these states are much too close. 

The sensitivity and time resolution of the micro- 
SQUID magnetometer |l^] allowed to study very small 
V15 crystals in good contact with the thermal bath. 
Down to 50 mK non-equilibrium behavior was never- 
theless observed at fast sweeping rate (up to dB /dt = 
0.7 T/s). Few hysteresis loops are represented Fig. ||a 
and Fig. ^a (only the positive part, the other one be- 
ing rigorously symmetrical). When the field increases, 



coming from negative values, the magnetization passes 
through the origin of the coordinates, reaches a plateau 
and then approaches the saturation. This leads to a 
winged hysteresis loop characterized by the absence of 
irreversibility near zero field. The wings depend sensi- 
tively on temperature and sweeping field rate. As an 
example, in Fig. |^ where three hysteresis loops are pre- 
sented (T=0.1 K, 0.15 K, 0.2 K) for a given sweeping rate 
(0.14 T/s), the plateau is higher and more pronounced 
at low temperature. The same tendency is observed at 
given temperature and faster sweeping rate (Fig. ^J) . At a 
given temperature, the equilibrium magnetization curve 
can be approximated by the median of the two branches 
of a low sweeping rate hysteresis loop (e.g. in Fig. [| for 
dBo/dt = 4.4 mT/s). When compared to the equilibrium 
curve, a given magnetization curve shows that (i) in nega- 
tive fields the system is colder than the bath, (ii) the spin 
system remains colder until the magnetization curve in- 
tersepts the equilibrium curve, (iii) after the intersept the 
spin temperature overpass the bath temperature, (iv) at 
sufficiently high fields (about 0.5 T) the system reaches 
the equilibrium. 

In order to interpret this magnetic behavior of the 
V15 molecules, we will analyse how the level occupa- 
tion numbers vary in this two level system (see Fig. ^a 
inset) when sweeping an external field. In the bare 
Landau-Zener (LZ) model |18 1?| the probability for the 
|l/2, -1/2) <-> |l/2, 1/2) transition is: 



exp 



(0.2) 



Taking the typical value r = 0.1 T/s and the zero-field 
splitting Ao = 0.05 K |20|pl| ], one gets a ground state 
switching probability very close to unity: the spin occu- 
pancy of the two level system is not modified when the 
field is reversed. In such an LZ transition, the plateaus 
of Fig. [I| should decrease if the sweeping rate increases, 
which is contrary to the experiments. Not surprizingly 
the 5=1/2 spin system is not isolated: LZ transitions 
are non-adiabatic due to spin couplings to the environ- 
ment. 

The mark of the V15 system is that this coupling is act- 
ing also near zero-field because Uuj « A is of the order of 
the bath temperature. The spin temperature T$ is such 
that ni/n2 — exp(A#/fcsT!s)> where Ah is the two lev- 
els field-dependent separation, and ni ; 2(^i,2eg) the out of 
equilibrium (equilibrium) level occupation numbers. In 
the magnetization curves at 0.1 K (Fig. ||a, ||a), the spin 
temperature is significantly lower than the bath temper- 
ature T (ni > ni eq , Ts < T) between —0.3 T (when 
the magnetization curve departs from the equilibrium 
one) and 0.15 T (the field at which the magnetization 
curve intersects the equilibrium one). After this inter- 
sept Ts is larger than the bath temperature (ni < ni eq , 
T s > T), and at sufficiently high fields (about 0.5 T) it 
reaches the equilibrium value (ni = n\ eqi T5 = T). Note 
that the magnetization curves measured between — 0.7 T 
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and 0.02 T at fast sweeping rates (0.07 and 0.14 T/s) 
are nearly the same, suggesting weak exchange with the 
bath, i.e. nearly adiabatic demagnetization. 

Indeed, as we will show below, the silver sample holder 
is not sufficient to maintain the phonon temperature 
equal to the temperature of the bath. This is because 
in V15 below 0.5 K, the heat capacity of the phonons 
C p h is very much smaller than that of the spins Cs, so 
that the energy exchanged between spins and phonons 
will very rapidly adjust the phonons temperature T p h 
to the spin one T s . Thus, we can see the spin sys- 
tem and the phonons as a single coupled system (quan- 
tum non-adiabatic LZ transitions) in weak exchange with 
the external bath (thermodynamical adiabatic demag- 
netization). The spins energy is transfered from the 
spins only to those phonons with Tilo — Ah (within 
the resonance line width). The number of such lattice 
modes being much smaller than the number of spins, 
energy transfer between the phonons and the sample 
holder must be very difficult, a phenomenon known as 
the phonon bottleneck |0. Following (HI, the num- 
ber of phonons per molecule available for such resonant 
transitions is tlt = J AuJ a(u>)duj / (e~xp{huj/kT) — 1), where 
a(ui)dui — 3Vu! 2 dcu/ (2tt 2 v 3 ) = number of phonon modes 
between lo and us + duj per molecule of volume V, v is 
the phonon velocity and Aus is the transition lincwidth 
due to fast hyperfine field fluctuations (they broden both 
energy levels). Taking the typical values v re 3000 m/s, 
T re 10" 1 K and Au re 10 2 MHz || we find n T of the 
order of re 10 -6 to 10 -8 phonons/molccule. Such a small 
number of phonons is very rapidly absorbed, burning a 
hole of width Aw in the spin-phonon density of states at 
the energy has — Ah 11 1 ■ If this spin-phonon density 
of states does not equilibrate fast enough, the hole must 
persist and move with the sweeping field, leading to a 
phonon bottleneck. 

4. Quantitative approach and numerical calcu- 
lations. 

Now this description will be made quantitative. For a 
given splitting Ah, the time evolution of the two levels 
populations ni,2 and of the phonon numbers riT ph at T p h 
obeys the set of two differential equations pi : 



-hi — h 2 — P\ 2 n\ ~ P 2 \n 2 

^T ph = -(n>T ph - n T )/T ph - Piim + P 2 in 2 , 



(0.3) 



where P\ 2 _ 2 i are the transition probabilities between 
the two levels (they are themselves linear functions of 
n T ph ) an d Tph re L/2v is the phonon-bath relaxation 
time (L is the sample size). Using the notations x — 
(n 1 -n 2 )/(n leq -n 2 e q ),y = (n Tph -n T )/(n T + n/2) with 
11 = Jaw CT M dw we S et: 



x = (1 — x — xy)jr\ 

v = -y/fph + bx, 



(0.4) 



where b = Cs/C p h and 1/tj = P\2 + P21 the direct spin- 
phonon relaxation time. By solving numerically this sys- 
tem for typical values, e.g. t\ = 10 -2 s, r p h < 10~ 6 s, 
b > 10 5 , we can see that T p h — > Tg ^ T (phonon bottle- 
neck) very rapidly, as expected. This leads to y = l/x — l 
and the second equation of the differential system be- 
comes x = (x — x 2 )/(l + bx 2 )/T p h. In the limit b >> 1 
(in our case b re 10 8 — 10 10 ) this equation has the solution: 



t/bTp h 



x + H(x - l)/(x - 1)), (0.5) 



where £0 = x(t = 0) and br p h is the spin-phonon relax- 
ation time (Tph = T$ —> T). When the system is not 
far from equilibrium (x ~ 1), we get an exponential de- 
cay of the magnetization, with the same time constant 
th = br p h- For a spin 1/2 system [^0[: 



th 



tanh 2 (A ff /2/c B T) 



A 2 

H 



(0.6) 



with a = 2-7T h v NTph/3Auj (N the molecule density). 

The dynamical magnetization curves calculated in this 
model are given Fig. |3|b and Fig. ^b. We started from 
equilibrium (xq = 1) in large negative fields. Then we 
let the system relax for a ver y short time St and we cal- 
culated x(8t) using Eq. 3.5. This value was taken as 
the initial value for the next field (the field step is r5t). 
The parameters have been chosen to mimic the measured 
curves of Fig. ||a and Fig. [|a |Q . The obtained similarity 
supports the possibility of the phonon bottleneck effect 
at the timescale of a few 0.1 s. In the Fig. [|a inset, we 
show the variation of the calculated spin-phonon temper- 
ature T s for T = 0.1 K and r = 0.14 T/s. After a cooling 
in negative fields, we can note a linear variation in the 
plateau region (small positive fields), where n\ln 2 re est.. 
The slope of this quasi-adiabatic linear region gives the 
plateau position and varies with the bath temperature 
and sweeping rate. In the Fig. ^b inset we show the calcu- 
lated field evolution of the number of phonons at energy 
%uj = Ah at equilibrium (T p h — Ts — T, dashed line) 
and out-of-equilibrium (riT ph = n-T=T SJ T — 0.14 T/s, 
black line) . The difference between the two curves (thick 
segment Auj) suggests the moving hole in the phonon dis- 
tribution, while their intersection gives the plateau inter- 
cept of the equilibrium magnetization curve (above which 
the hole dissapears and T p h = Ts > T). Let note that in 
zero field the system is out-of-equilibrium even if magne- 
tization passes through the origin of coordinates (without 
a barrier, the switch between +1/2 and —1/2 follows the 
level structure shown Fig. || inset ). At larger fields, in 
the plateau region, ni/n 2 re est. at timescales shorter 
than th = br p h (Eq. |0.6[ ), even after the plateau crosses 
the equilibrium curve. Equilibrium is reached when th 
becomes small enough. 

Furthermore, we measured the relaxation of the mag- 
netization of our crystal at different fields and temper- 
atures, along the plateau region. The relaxation curves 
compared well to exponential decay and the obtained re- 
laxation times are presented in Fig. |^a. The comparison 
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with those calculated (Fig. 
noted that a direct fit to Eq. 
values for a and A (« 0.4- 

5. Conclusions 



pp) is acceptable. But we 
|0.5| would necessitate larger 
0.6 sK 2 and w 0.2- 0.3 K). 



In conclusion, this dynamical study of a single crys- 
tal of non-interacting V15 molecules with spin 1/2 is an 
example of a non-adiabatic Landau-Zener model with 
dissipation. Due to relatively large splitting A at the 
5=1/2 anticrossing in V15 molecules (i) the probability 
for ground-state adiabatic transitions is nearly equal to 
one and (ii) the spin system absorbes phonons during the 
Landau-Zener transition creating a hole in their distribu- 
tion. The time and field evolution of this hole generates 
a "butterfly" hysteresis loop quite different from the one 
of high spin molecular magnets with large barrier and in- 
finitesimal Ao/hsT ratio (no phonons are available at the 
anticrossing). The effects presented in this paper seem 
to be a mark of molecules with low spin. 
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FIG. 1. The geometry and the spin coupling scheme in V15 
complex. There are five antiferromagnetic exchange constants 
leading to a ground state S — 1/2 in fields below 2.8 T and 
S — 3/2 for greater fields (here such a field is pointing down). 




FIG. 2. The agreement between the calculated curve at 
thermodynamical equilibrium (line) and the experimental 
magnetization data is quite good for an isotropic exchange 
Jo = -2.445 K. For J x = Jy = -2.75 K, J z = -2.2 K 
and the field applied in the triangle plane (dashed line) one 
observes that the calculated second flip is slightly displaced 
but its slope is still greater than the experimental one. In 
inset are represented the energy levels vs. applied field in the 
effective triangle approach. When the applied field reaches 
—3Jo/{2g/j,B) the ground state switch from S = 1/2 to 
S = 3/2. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
B (T) 

FIG. 3. Measured (a— top) and calculated (b— bottom) hys- 
teresis loops for three temperatures and for a given field 
sweeping rate 0.14 T/s. The plateau is more pronounced at 
low T. The inset is a schematic representation of a two-level 
system Sz = ±1/2 with repulsion due to non-diagonal matrix 
elements. In a swept field the switching probability P is given 
by the Landau-Zener formulae (see text). The two levels are 
broadened by the hyperfine fields and the absorption or the 
emission of phonons can switch the polarization state of spins. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
B (T) 

FIG. 4. Measured (a— top) and calculated (b— bottom) hys- 
teresis loops for three field sweeping rates at T = 0.1 K. The 
observed plateau is more pronounced at high sweeping rate. 
The equilibrium curve can be approximated by the median 
of the two branches of the low sweeping rate hysteresis loop 
(dotted curve). In the top inset is plotted the spin and phonon 
temperature T s = T ph for T = 0.1 K and r = 0.14 T/s, when 
the field is swept from negative values. Ts decreases until 
zero-field and then increases linearly within the plateau re- 
gion. Then it overpasses the bath temperature to finally reach 
the equilibrium. In the bottom inset the calculated number of 
phonons with hu = Ah is plotted vs. the sweeping field mod- 
ulus (note the arrows) at equilibrium (T p h — Ts — T, dashed 
line) and out-of-equilibrium (tit h = ut=t s , r = 0.14 T/s, 
black line) . The difference between the two curves (thick seg- 
ment Acj) suggests the moving hole in the phonon distribu- 
tion, while their intersection gives the plateau intercept of the 
equilibrium magnetization curve. 
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FIG. 5. The relaxation times th, measured (a— top) and 
calculated (b— bottom, same parameters as in Fig. ^, ^b). 
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